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Abstract. The mid-Pliocene Warm Period (mPWP) bears 
many similarities to aspects of future global warming as pro- 
jected by the Intergovernmental Panel on Climate Change 
(IPCC, 2007). Both marine and terrestrial data point to high- 
latitude temperature amplification, including large decreases 
in sea ice and land ice, as well as expansion of warmer 
climate biomes into higher latitudes. Here we present our 
most recent simulations of the mid-Pliocene climate using 
the CMIP5 version of the NASA/GISS Earth System Model 
(ModelE2-R). We describe the substantial impact associated 
with a recent correction made in the implementation of the 
Gent-Mc Williams ocean mixing scheme (GM), which has a 
large effect on the simulation of ocean surface temperatures, 
particularly in the North Atlantic Ocean. The effect of this 
correction on the Pliocene climate results would not have 
been easily determined from examining its impact on the 
preindustrial runs alone, a useful demonstration of how the 
consequences of code improvements as seen in modem cli- 
mate control mns do not necessarily portend the impacts in 
extreme climates. 

Both the GM-corrected and GM-uncorrected simulations 
were contributed to the Pliocene Model Intercomparison 
Project (PlioMIP) Experiment 2. Many findings presented 
here corroborate results from other PlioMIP multi-model 
ensemble papers, but we also emphasise features in the 
ModelE2-R simulations that are unlike the ensemble means. 
The corrected version yields results that more closely resem- 
ble the ocean core data as well as the PRISM3D reconstruc- 
tions of the mid-Pliocene, especially the dramatic warming 
in the North Atlantic and Greenland-Iceland-Norwegian Sea, 


which in the new simulation appears to be far more realis- 
tic than previously found with older versions of the GISS 
model. Our belief is that continued development of key phys- 
ical routines in the atmospheric model, along with higher res- 
olution and recent corrections to mixing parameterisations in 
the ocean model, have led to an Earth System Model that will 
produce more accurate projections of future climate. 


1 Introduction 

Climate simulations with versions of the NASA/GISS gen- 
eral circulation models have been used to explore the 
Pliocene as a potential future climate analogue since NASA 
and the USGS partnered on a data-model comparison project 
in the early 1990s (Rind and Chandler, 1991; Chandler and 
Rind, 1992; Chandler et al., 1994; Poore and Chandler, 
1994). However, as with numerous other climate modelling 

studies, the GISS model has often underestimated the high 
degree of polar amplification seen in Pliocene proxy data, 
particularly in the North Atlantic Ocean, without resorting 
to levels of CO 2 that are extreme and which are not sup- 
ported by other studies (Hansen and Sato, 2012; Pagani et 
al., 2010; Seki et al., 2010). To further complicate the issue, 
higher CO 2 levels lead to low-latitude sea surface tempera- 
ture increases that are not supported by an increasing num- 
ber of ocean core analyses from tropical locations, a find- 
ing which largely validates prior assertions that the Pliocene 
tropical sea surface temperatures (SSTs) were little changed 
from the modern. Still, the ability to compare warm-climate 
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simulations with data-supported global reconstructions is an 
advantage that is not afforded by future climate change sce- 
narios. Many groups, therefore, continue to pursue a bet- 
ter understanding of the middle Pliocene as it provides cli- 
mate scientists with one of the few warm-climate scenarios 
in which global and high-latitude temperatures were as warm 
as IPCC’s future climate projections, and has continental and 
ocean basin configurations that approximate modern geogra- 
phy. 

The experiments examined in this study were performed 
in conjunction with the Pliocene Model Intercomparison 
Project (PlioMIP) Experiment 2 (Phase 1), the coupled 
ocean-atmosphere simulation. A description of the experi- 
ment design can be found in Haywood et al. (2010, 2011) and 
the description of the large-scale features of the multi-model 
ensemble is found in Haywood et al. (2013). Specific model 
characteristics that affect experiment design and which are 
unique to the NASA/GISS modelling effort are described be- 
low. 


2 PlioMIP experiment design 

For PlioMIP, a set of environmental reconstructions of the 
mid-Piacenzian Stage was compiled to form the PRISM3D 
reconstruction (Dowsett et al., 2010a). These reconstructions 
are distributed by the US Geological Survey as a series of 
uniformly gridded (2° x 2°) datasets, and constitute the pre- 
ferred PlioMIP experiment protocol (Haywood et al., 2010, 
2011). The adaptation of these datasets into boundary condi- 
tions for use with GISS ModelE2-R, as well as the creation 
of required related inputs, is described below. 

2.1 Land/ice topography and ocean bathymetry 

A reconstructed topography and land/sea mask for the 
Pliocene (Sohl et al., 2009), reflecting a 25-m rise in sea 
level compared to the modern, was developed as part of the 
PRISM3D effort. Elevation and areal distribution informa- 
tion for static representations of the Greenland Ice Sheet and 
East Antarctic Sheet (Hill et al., 2007) are incorporated into 
the topography, as the ice sheets must essentially be repre- 
sented as “big white plateaus” for models that do not have dy- 
namic land ice capabilities. The original 2° x 2° PRISM3D 
gridded topographic dataset was regridded in two ways for 
use with GISS ModelE2-R: once to 2° latitude x 2.5° longi- 
tude for the atmosphere component of the model, and further 
to 1° latitude x 1.25° longitude for the ocean component; for 
the land/sea mask, both atmosphere and ocean grids are non- 
fractional. The topographic relief from the resulting higher- 
resolution dataset is used directly as the input boundary con- 
dition, and is not applied as an anomaly relative to the mod- 
ern as described in Haywood et al. (2010). Corrections to 
the resulting land-sea masks at both resolutions were made 
by hand to ensure that continental outlines were consistent 


after regridding, that larger islands did not disappear as a re- 
sult of the interpolation used in the regridding process, and 
that narrow ocean passages that existed in the Pliocene re- 
mained open. Note that in ModelE2-R, a straits parameteri- 
sation is used to maintain ocean flow through grid cell loca- 
tions where straits cannot be resolved at the 1° x 1.25° res- 
olution of the land/sea mask in the ocean component of the 
model. The changes to the land/sea mask introduced by the 
25-m sea level rise removed the need to parameterise certain 
straits, since the straits could now be resolved; in these loca- 
tions, Pliocene ocean waters flow freely. (See Table 1 for a 
comparison to the modern). 

Note that the entire West Antarctic Ice Sheet is absent in 
our reconstruction (as per Pollard and DeConto, 2009), creat- 
ing an “Ellsworth Passage” that separates the modern Antarc- 
tic Peninsula from the rest of Antarctica. The Pliocene ocean 
bathymetry in this region is uncertain, since we do not know 
what the thickness of any previous West Antarctic Ice Sheet 
may have been, nor whether the region was still undergoing 
glacio-isostatic adjustment after the loss of such ice. The cur- 
rent depth to bedrock in the region varies from 500 to 2000 m 
(Lythe et al., 2000); we used a uniform depth of approxi- 
mately 500 m as a reasonable estimate for maximum ocean 
depth through the passage during the Pliocene. 

Ocean bathymetry for the Pliocene was otherwise adapted 
from the modern, with most modifications made to accom- 
modate continental edges flooded by the increase in sea level. 

2.2 Riverflow and continental drainage 

The river drainage system of the continents for the Pliocene 
experiments is similar to that of modem geography, but ex- 
ceptions exist that cannot be prescribed based on the avail- 
able paleogeographic evidence. For example, adjustments re- 
lated to the coastal slope changes are required due to the sea 
level rise of 25 m, which was prescribed as consistent with 
reductions of ice sheet volume. Using the topographic ele- 
vation boundary condition array, and working inward from 
continental edges, we calculate the slope of each continental 
grid cell in eight directions (four sides, four corners). Runoff 
is then removed from each cell in the direction of maximum 
slope, tracing a route back to the coast. For coastal grid cells 
that have more than one border adjacent to an ocean grid cell, 
mnoff crosses the coastal grid cell on the same trajectory as 
in the adjacent inland grid cell. 

Glacial ice melts directly from the Greenland and East 
Antarctic ice sheets, and enters the surface ocean in pre- 
scribed cells wherever the edges of the ice sheets coincide 
with continental edges. 

2.3 Ocean temperatures and salinity 

Both the sea surface and deep ocean temperature datasets 
provided with PRISM3D (Dowsett et al., 2009, 2010b) 
were regridded to 1° latitude x 1.25° longitude horizontal 
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Table 1. Use of straits parameterisations in GISS ModelE2-R - Modern vs. Pliocene. 


Strait Name 

Geographic Location 

Depth (m) 

Width (m) 

Modern 

Pliocene 

Dolphin & Union 

Between Victoria Island and mainland Canada, 
Northwest Territories/Nunavut 

56 

32 000 

parameterised 

explicit 

Dease 

Between Victoria Island and the Kent Penin- 
sula, Nunavut, Canada 

56 

25 000 

parameterised 

explicit 

Fury & Hecla 

Between Baffin Island and the Melville Penin- 
sula, Nunavut, Canada 

30 

15 000 

parameterised 

parameterised 

Nares 

Between Ellesmere Island, Canada and Green- 
land 

202 

30 000 

parameterised 

parameterised* 

Gibraltar 

Between Spain and Morocco 

280 

15 000 

parameterised 

explicit 

English 

Between England and France 

30 

35 000 

parameterised 

explicit 

Dardanelles 

Northwestern Turkey 

30 

1000 

parameterised 

parameterised 

Bosporous 

Northwestern Turkey 

30 

2000 

parameterised 

explicit 

Bab al-Mandab 

Connects Red Sea and Gulf of Aden 

140 

25 000 

parameterised 

parameterised 

Malacca 

Between Malay Peninsula and Sumatra, In- 
donesia 

30 

40 000 

parameterised 

parameterised 

Selat Sunda 

Between Java and Sumatra, Indonesia 

30 

25 000 

parameterised 

parameterised 


* The Nares Strait is divided into two parts to accommodate Pliocene geography. 


resolution; the vertical resolution of the deep ocean temper- 
ature dataset was maintained at 33 layers. Wherever there 
were missing values within the regridded temperature dataset 
with respect to our new ocean bathymetry, we used a “nearest 
neighbour” approach to fill the gaps, always maintaining the 
integrity of the vertical temperature profile. Salinity values 
for the Pliocene were derived from modern salinity values 
(Conkright et al., 1998), with missing values filled in in the 
same fashion as the temperatures. 

2.4 Biome distribution mapping to GISS ModelE2-R 
vegetation categories 

GISS ModelE2-R uses a vegetation scheme that includes 
eight vegetation types, plus land ice as an additional type. 
Each type is used in to define or modify certain physical 
characteristics of the land surface and ground, such as vis- 
ible and near-infrared surface albedo, the water holding ca- 
pacity of soil layers, transpiration rates and snow masking 
depths. The GISS vegetation categories were originally dis- 
tilled from a more detailed global vegetation compilation of 
Matthews (1983, 1984), and subsequently have been modi- 
fied to better reflect agricultural coverage and its impact. Ir- 
rigation effects are not included. The GISS model’s radiation 
code also accounts separately for subgrid-scale fractions of 
bare ground in each cell. 

The prescribed Pliocene biome distributions for PlioMIP 
were developed by the British Antarctic Survey (BAS biomes 
and mega-biomes; Salzmann et al., 2008), with each partic- 
ipating modelling group being responsible for developing a 
method of translation appropriate for their model’s vegeta- 
tion parameterisation. To make this translation for the GISS 
model, we regridded the BAS modem mega-biome distribu- 
tion to 2° latitude x 2.5° longitude (the standard resolution 


Table 2. Correlation of BAS mega-biomes with GISS ModelE2-R 
modern vegetation classes. 


BAS Mega-Biome Type 

GISS Vegetation Type 

Tropical Forest 

Rainforest 

Warm-Temperate Forest 

Deciduous forest 

Savanna/dry woodland 

Grassland 

Grassland/dry shrub 

Grassland 

Desert 

Desert 

Temperate forest 

Deciduous forest 

Boreal forest 

Evergreen forest 

Tundra 

Tundra 

Dry tundra 

Tundra 


for the ModelE2-R vegetation input), and then compared 
the BAS distribution cell by cell against the GISS modem 
vegetation distribution used for ModelE2-R. We then used 
a frequency distribution chart plotting BAS vs. GISS types 
to identify dominant links between types (see Table 2), and 
the BAS Pliocene mega-biome distribution was translated ac- 
cordingly, with non-fractional values assigned to each cell. 
Note that under this translation scheme, the GISS vegetation 
types shmb/grassland and tree/grassland have no clear corol- 
lary to the BAS mega-biomes, and instead are represented by 
grassland alone. 

2.5 NASA/GISS ModelE2-R 

The climate model used for these simulations was developed 
at the NASA Goddard Institute for Space Studies (GISS) and 
is the CMIP5/PMIP3 version that will be archived as GISS 
ModelE2-R for IPCC AR5 and PMIP. The suffix (-R) is used 
by GISS to identify the specific ocean version in the coupled 
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model, in this case the Russell ocean model. The physics 
and parameterisations of the model are described in Schmidt 
et al. (2006) and updates for AR5 are being described in 
Schmidt et al. (2012, 2013). The most recent documenta- 
tion is also available through the NASA/GISS website. This 
model is part of a multi-decade code lineage developed at 
NASA/GISS and which has been referred to in the literature 
over the years as Model IT, si 1995, si 1997 and si2000 (e.g., 
Hansen et al., 2000, 2002) and afterward as ModelE. Hence- 
forth, in this paper we will refer to the model as ModelE2-R, 
or ModelE2 if we are describing the atmosphere-only ver- 
sion. 

ModelE2-R calculates temperature, pressure, winds and 
specific humidity as prognostic variables, using the conser- 
vation equations for mass, energy, momentum and moisture. 
The standard configuration produces global climate simula- 
tions at a latitude x longitude resolution of 2.0° x 2.5° in the 
atmosphere, 1.0° x 1.25° in the ocean, and includes 40 lay- 
ers in the atmosphere and 32 layers in the fully coupled dy- 
namic ocean (Russell et al., 1995). ModelE2 (atmosphere 
only) uses second-order differencing schemes in the momen- 
tum and mass equations and a quadratic-upstream scheme for 
heat and moisture advection, which implicitly enhances an 
absolute model resolution of 2.0° x 2.5° grid to 0.7° x 0.8° 
(Schmidt et al., 2006, Table 1). The radiation physics in- 
cludes calculations for trace gas constituents (CO 2 , CH 4 , 
N 2 O, CFCs, O 3 ) and aerosols (natural and anthropogenic) 
and is capable of simulating the effects of large forcing 
changes in constituents such as volcanic aerosols and green- 
house gases. The forcings are assigned at startup or can be 
altered transiently throughout an experiment. In addition, a 
parameterised gravity wave drag formulation incorporates 
gravity- wave momentum fluxes that result from flow over to- 
pography and deformation in baroclinic systems. The gen- 
eration, propagation and drag are all a function of the calcu- 
lated variables at each grid box for the various vertical levels. 
New to this model is the capability of running with interac- 
tive atmospheric chemistry, aerosols and dust concentrations, 
as well as the incorporation of a model-calculated (rather 
than parameterised) first aerosol indirect effect. ModelE2-R 
also employs a prognostic cloud water parameterisation (Del 
Genio et al., 1996; Schmidt et al., 2006, 2012) that represents 
all-important microphysical processes. 

2.6 Corrections to the ocean mixing scheme 

It is important to emphasise that the simulations described 
here as “GM-CORR” utilise a correction to the Gent- 
McWilliams parameterisation in the ocean component of 
the coupled climate model. In prior implementations of the 
mesoscale mixing parameterisation in GISS ModelE, which 
like many ocean models uses a unified Redi/GM scheme 
(Redi, 1982; Gent and McWilliams, 1990; Gent et al., 1995; 
Visbeck et al., 1997), a miscalculation in the isopycnal 
slopes led to spurious heat fluxes across the neutral surfaces, 


resulting in an ocean interior generally too warm, but with 
southern high latitudes that were too cold. A correction to re- 
solve the problem was made for this study, and it will also 
be employed in all subsequent versions of ModelE2-R going 
forward. In applying the correction the new code also uses 
a minimum eddy diffusivity of 600 m 2 s -1 in the mesoscale 
eddy formulation. Since some multi-model ensembles have 
included the previous version we also present here a compar- 
ison of some of the more dramatic differences between the 
old and new simulations. In this paper, we refer to the uncor- 
rected version of the model as “GM-UNCOR”. No other dif- 
ferences exist between the two simulations beyond the cor- 
rections to the mixing scheme. All boundary conditions, ini- 
tial conditions and model physics are identical and there were 
no changes made to the land-sea masks or subgrid-scale strait 
definitions. 


3 Simulation results 

All results presented in the following sections are climatolog- 
ical averages for years 921-950 of our Pliocene and Prein- 
dustrial simulations. The sea surface temperatures used in 
the Preindustrial control run are from 1876-1885. As noted 
above in Sect. 2.6, the main results from the GISS ModelE2- 
R now include a correction to the ocean mixing scheme. 
Pliocene and Preindustrial runs including the correction are 
labeled “GM_CORR”; results that are from the uncorrected 
version of the model are labeled “GMJJNCOR”. 

3.1 Sea surface and surface air temperature response 

The most fundamental mismatch between warm climate pa- 
leoclimate simulations and paleoclimate proxy data has been 
the inability of simulations to achieve both a reduced pole- 
to-equator temperature gradient along with acceptable mag- 
nitudes of polar and tropical temperature amplification. Ex- 
cessive greenhouse gas (GHG) levels have traditionally been 
required for climate models to achieve polar temperature am- 
plification that is in the vicinity of the observed levels for the 
warmer time periods in Earth history. However, that same 
GHG excess tends to cause tropical temperatures that are un- 
realistically high. 

Pliocene simulations are no different in this respect, in- 
cluding those for PlioMIP Experiment 2 (Haywood et al., 
2013). However, PlioMIP Pliocene experiments prescribe the 
use of 405ppmv CO 2 , with other greenhouse gases held at 
preindustrial levels. Though 405 ppmv is actually on the up- 
per end of proxy-based estimates for the Pliocene (Kiirschner 
et al., 1996; Raymo et al., 1996; Pagani et al., 2010), that 
level is relatively low from the perspective of other warm Ter- 
tiary and Mesozoic climates. The problem of tropical tem- 
perature amplification in Pliocene simulations is thus less 
of a conundrum than for most past warm periods. Given 
that there are still large uncertainties associated with pre- 
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ice core CO 2 estimates, and given that the tropical energy 
balance is very sensitive to even small CO 2 changes, the 
simulated tropical SSTs are within the envelope of uncer- 
tainties. All the same, the warming in mid- to high lat- 
itudes is supplemented substantially by the direct forcing 
and feedbacks associated with the additional carbon dioxide; 
405 ppmv is, after all, 45 % higher than the 280 ppmv value 
used in most preindustrial control runs. Of course, high- 
latitude temperatures are impacted by the specified sea ice 
and ice sheet distributions, which are substantially decreased 
compared with the present day (Hill et al., 2007; Dowsett et 
al., 2010b); however, compared to the CMIP3 model results 
and previous GISS model simulations of the Pliocene (e.g., 
Chandler et al., 1994; Shukla et al., 2009, 2011; Haywood 
et al., 2013), the GM.CORR version of the GISS ModelE2- 
R shows compelling improvement. In the final analysis, the 
corrected simulation results for PlioMIP Experiment 2 yield 
sea surface temperatures in the North Atlantic that are closer 
to PRISM3D core data than those that we previously con- 
tributed to the multi-model ensemble of Dowsett et al. (2012) 
(see cool bias in the North Atlantic region of Fig. 3e from 
Dowsett et al., 2012). 

3.2 Global and zonal average temperatures 

The global warming compared to preindustrial simula- 
tions is similar in both Pliocene simulations (Table 3) at 
2.25 °C for the Pliocene GM.CORR and 2. 12 °C for Pliocene 
GM_UNCOR. Zonal average sea surface and surface air tem- 
perature anomalies (SST and SAT in Fig. 1) for the Pliocene 
GM.CORR simulation peak over the regions of sea ice loss 
and ice sheet reduction. The surface air temperatures are 
slightly warmer than the multi-model ensemble averages 
for all PlioMIP models in both the Northern (+8.2 °C) and 
Southern (+9.1 °C) Hemispheres (see Haywood et al., 2013). 
At low latitudes the character of the SAT temperature anoma- 
lies are similar between the GM_CORR and GM_UNCOR 
versions of the model. However, the ocean mixing correction 
has decreased zonal average temperatures in the extratropics 
of the Southern Hemisphere while increasing temperatures 
in the Northern Hemisphere. The significant exception is in 
the northern high latitudes, where surface air temperatures in 
GM_CORR are warmer than in GM_UNCOR by at least 2 °C 
at all latitudes north of 50° N. 

High latitude amplification of the sea surface temperature 
anomaly is asymmetrical, with more than twice as much 
warming in the Northern Hemisphere (+4.5 °C) as in the 
Southern Hemisphere (+2.2 °C) in GM_CORR. The peak of 
the zonal average SST falls at 65° N, at the high end of the 
range of PlioMIP models, and slightly exceeds the PRISM3D 
data for that latitude. Sea surface temperatures in GM_CORR 
are up to 9 °C warmer in the North Atlantic than in the un- 
corrected model at their peak (Fig. 2); in fact, zonal aver- 
age SSTs in the corrected model are warmer than at any 
other location. In contrast, at similar latitudes the peak SST 
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Fig. 1. Zonal average anomalies of surface air temperature (in blue) 
and sea surface temperature (in red). The simulated high-latitude 
amplification of surface air temperature and the peak SSTs around 
65° N are more similar to proxy data than those produced by the 
earlier GMTJNCOR version of the GISS ModelE2-R, but tropical 
temperature anomalies in both versions are still slightly higher than 
data interpretations for the mid-Pliocene. 


anomalies in the uncorrected model actually showed cool- 
ing compared to their preindustrial counterpart. The regional 
aspects of this large difference caused by the ocean mixing 
correction will be discussed below. 

The correction to the ocean mixing scheme clearly is 
a substantial improvement on the previous GM_UNCOR 
Pliocene simulation, which was too cool, zonally averaged, 
at high northern latitudes due to the large cooling in the North 
Atlantic. The peak SST warming in the Southern Ocean 
and the 1-2 °C zonal average increase across southern mid- 
latitudes are broadly similar to proxy data in the region, 
which show little longitudinal variation and yield tempera- 
ture increases in the same range. 

Tropical zonal average SAT anomalies are approximately 
1.1 °C and sea surface temperatures in the tropics are simi- 
larly amplified, with a zonal average warming of 1.2 °C just 
north of the equator. Both corrected and uncorrected models 
yield broadly similar results in the tropics. The values fall on 
the low side of the range of PlioMIP model results for the 
region, but remain a bit warmer than the PRISM3D recon- 
struction, which is barely elevated above late 20th century 
tropical temperatures within 10° of the equator. 

3.3 Regional temperature changes: North Atlantic and 
Arctic oceans 

Using a four-model subset of the PlioMIP group of cli- 
mate models (including the Pliocene GM_UNCOR version 
of ModelE2-R), Dowsett et al. (2012) generated a mean SST 
field to compare with the PRISM3 ocean core proxy data. 
They reported several regional differences between proxy 
data and simulation results, but the most pronounced mis- 
match of models and data was in the underestimate of North 
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Table 3. Global mean climate variables for Pliocene and Preindustrial simulations 



GM_CORR 

GMJJNCOR Difference 

GM_CORR 

GMJJNCOR Difference 


Surface Air Temperature (°C) 


H 2 O of Atmosphere (mm) 


Pliocene 

16.21 

16.29 

-0.08 

25.99 

25.91 

0.08 

Preindustrial 

13.96 

14.17 

-0.21 

22.74 

22.66 

0.08 

Difference 

2.25 

2.12 

0.13 

3.25 

3.25 

0.00 


Planetary Albedo (%) 


Precipitation (mm day 1 ) 


Pliocene 

28.87 

28.72 

0.15 

3.34 

3.34 

0.00 

Preindustrial 

29.71 

28.55 

1.16 

3.17 

3.17 

0.00 

Difference 

0.84 

0.17 

-1.01 

0.17 

0.17 

0.00 


Snow Coverage (%) 


Evaporation (mm day 1 ) 


Pliocene 

7.67 

7.3 

0.37 

3.34 

3.34 

0.00 

Preindustrial 

9.96 

9.26 

0.70 

3.17 

3.17 

0.00 

Difference 

-2.29 

-1.96 

-0.33 

0.17 

0.17 

0.00 


Snow Depth (mm^O) 


Total Cloud Cover (%) 


Pliocene 

12.8 

12.12 

0.68 

59.51 

59.28 

0.23 

Preindustrial 

10.35 

9.72 

0.63 

61.76 

61.51 

0.25 

Difference 

2.45 

2.40 

0.05 

-2.25 

-2.23 

-0.02 


Global Ocean Ice Cover (%) 


Moist Convective Cloud Cover (%) 

Pliocene 

3.14 

2.72 

0.42 

4.2 

4.25 

-0.05 

Preindustrial 

4.71 

3.96 

0.75 

4.08 

4.15 

-0.07 

Difference 

-1.57 

-1.24 

-0.33 

0.12 

0.10 

0.02 


Arctic Ice Cover (%) 


Southern Ocean Ice Cover (%) 


Pliocene 

3.40 

4.04 

-0.64 

2.87 

1.41 

1.46 

Preindustrial 

5.71 

5.29 

0.42 

3.72 

2.62 

1.1 

Difference 

-2.31 

-1.25 

-1.06 

-0.85 

-1.21 

0.36 


Atlantic temperature change by the models (see Fig. 3e, 
Dowsett et ah, 2012). The North Atlantic SST anomalies 
produced by the GISS model before applying the ocean 
mixing correction (GMJJNCOR) were by far the least rep- 
resentative of the proxy data, actually exhibiting cooling, 
whereas most models produce some warming in the North 
Atlantic (though not enough to match data). The results us- 
ing ModelE2-R with the corrected ocean mesoscale mixing 
scheme (GM.CORR) are much warmer in the North Atlantic 
(Fig. 2), actually exhibiting more warming in that region than 
any other model in the PlioMIP multi-model ensemble. Yet, 
GM.CORR still exhibits small regional cooling in the up- 
per mid-latitude Atlantic, a bias also seen in the preindustrial 
simulations as compared to observed SST climatologies of 
the late 19th century (see Dowsett et al., 2012 Fig. 3d). 

Despite minor differences in the zonal average anomalies 
compared to previous GISS model simulations, the zonal val- 
ues are soundly within the range defined by the PlioMIP 
model ensemble for zonal temperature results. However, ex- 
amination of the regional temperature fields shows signifi- 
cant basin-to-basin differences and the substantial improve- 
ments that the ocean mixing scheme correction has yielded. 


In the previous GMJJNCOR version of the model (Fig. 3a) 
the dominant feature in the SST anomaly field was the strong 
cooling in the North Atlantic Ocean. The cooling peaked 
in the central North Atlantic around 60° N, but was broadly 
distributed throughout the North Atlantic, from the coast of 
Nova Scotia, to the Norwegian Sea, extending along the re- 
turn gyre all the way to Africa and into the tropical flow. The 
Arctic experienced a broad 1-2 °C warming, less than most 
other regions of the ocean. 

In stark contrast, the ocean mixing scheme correction 
leads to an intense SST warming in the Pliocene GM_CORR 
version of the model (Fig. 3b). The warm anomaly peaks 
in the Norwegian and Iceland sectors of the GIN Sea at 
+9.7 °C and extends uninterrupted to the south of Greenland 
into the Labrador Sea, where another peak exceeds +6 °C 
in the Atlantic Ocean east of Newfoundland. These very 
high anomalies do not, however, extend far into the Arctic 
Ocean, where the 1-2 °C warming is similar to that seen in 
the GMJJNCOR version of the model. The warm anomaly in 
the GIN Sea spreads southward, following the return flow of 
the North Atlantic gyre as it merges with the Canary Current, 
though the values are only about half as large in this region 
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SEA SURFACE TEMPERATURE 

Pliocene (GM_CORR) - Pliocene (GMJJNCOR) 
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Fig. 2. Corrections to the ocean mesoscale mixing scheme 
(GM_CORR) yield substantially warmer sea surface temperatures 
in the North Atlantic, closer to what ocean core data suggest. Al- 
though even GM_CORR still exhibits small regional cooling in the 
upper mid-latitude Atlantic, a bias also seen in the multi-model en- 
semble, preindustrial simulations as compared to observed prein- 
dustrial SST climatologies (1870-1900; see Dowsett et al., 2012 
Fig. 3d). 


as in the PRISM3D data (3^1 °C as opposed to 6-8 °C), 
whereas further north the peak values more closely resem- 
ble the ocean core data. In comparing the simulated warming 
to the core data poleward of 75° N, the model warming is 
also 2-4 °C too cool, but again the new simulation compares 
far more favourably to data reconstructions than the Pliocene 
simulation that was run with the uncorrected CMIP5 ver- 
sion of the model and is even an improvement compared to 
coarser resolution CMIP3 GISS ModelE-R simulations (not 
shown). Noticeably, the most pronounced negative temper- 
ature anomaly in the Pliocene ocean also lies in the North 
Atlantic in GM_CORR, though the geographic character is 
unlike the cool anomaly in the GMJJNCOR results (com- 
pare Fig. 3a and b). Figure 3b shows a region of cooling that 
is narrow in latitudinal extent, but spanning nearly the width 
of the Atlantic Ocean, tracing the regional expression and 
extent of the modem Gulf Stream Current. 

Delayed warming is characteristic of the North Atlantic 
region in many warm-climate coupled model simulations, 


SEA SURFACE TEMPERATURE 
Pliocene - Preindustrial (gmjjncor) 



Pliocene - Preindustrial (gm_corr) 
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Fig. 3. The sea surface temperature anomaly field from the Pliocene 
simulations compared to their respective preindustrial control runs. 
The anomalies show a generally warmer planet with high-latitude 
amplification of SSTs. The most apparent feature is the large tem- 
perature change in the North Atlantic and, especially, the dramatic 
warming that results when the corrected ocean mixing scheme is ap- 
plied. The North Atlantic and Southern Ocean warm anomalies are 
key fingerprints of the mid-Pliocene Warm Period as reconstructed 
by the US Geological Survey PRISM3D Project, suggesting that the 
corrected version of ModelE2-R (GM_CORR in the text) is a great 
improvement over the uncorrected version (GMJJNCOR) and is 
also an advance compared to simulations that used the coarser-grid 
CMIP3 version of GISS ModelE. 

including the PlioMIP ensemble simulations (compare 
Figs. 2a and 3a from Haywood et al., 2013) and the fu- 
ture climate change simulations documented in the IPCC’s 
Fourth Assessment Report (AR4; IPCC, 2007, see Fig. 10.8). 
In addition, the North Atlantic is one of the few regions 
where inter-model standard deviation exceeds the mean of 
the multi-model ensembles (Haywood et al., 2013; Meehl 
et al., 2007, see Figs. 10.8 and 10.9) revealing just how 
much our best models differ, and how great the uncertainty 
is regarding climate change in this critical region. It is not 
surprising then that simulations of the mid-Pliocene con- 
ducted using coupled ocean-atmosphere models have rarely 
achieved adequate simulations of the North Atlantic Ocean. 
Ocean core data and, therefore, PRISM3D reconstructions, 
have shown with considerable confidence that the North 
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Atlantic and northward extensions into the GIN Sea expe- 
rienced strong warming in the Pliocene - greater than any- 
where else in the global oceans and greater than in other ex- 
periments in the multi-model PlioMIP ensemble (Haywood 
et ah, 2013; see Fig. 3). Dowsett et al. (2012) point out that 
even accounting for a cool bias in the models compared to 
their own preindustrial control simulations, coupled models 
underestimate North Atlantic warming in the mid-Pliocene 
by 2-8 °C. 

3.4 Regional ocean temperature changes: North Pacific, 
Southern Ocean, and the Tropics 

Temperature change throughout the rest of the global oceans 
is muted by comparison with the North Atlantic and Arctic, 
but there are still significant anomalies that can be compared 
and contrasted to data reconstructions in other key oceanic 
regions. The North Pacific and Southern Ocean show warm- 
ing in the Pliocene in both the GM_CORR and GM_UNCOR 
simulations (Fig. 3a and b). Unlike the Atlantic anomalies, 
the sea surface temperatures in the Pacific and Southern 
Oceans are broad and peak below + 5 °C. In addition, the 
warming was 2^1 °C less once the ocean mixing scheme 
was corrected. More extreme warm anomalies near the Asian 
shoreline in the northwest Pacific are supported by ocean 
core data that Dowsett et al. (2012) give a confidence rat- 
ing of “Very High”. It is noteworthy that the temperature 
anomaly lessens poleward of 45° N latitude and into the 
northeastern Pacific Ocean in the GM_CORR model. Data 
from the Gulf of Alaska do not dispute that trend. Cores 
taken off the west coast of North America are given a “high- 
confidence” rating in PRISM project analyses and indicate 
levels of Pliocene warming that are not simulated by the cor- 
rected version of the model. 

While mid- to high latitude warming is evident in ocean 
cores throughout the Southern Ocean there is little longitu- 
dinal variation in the warming - despite the absence of the 
West Antarctic Ice Sheet in the Pliocene boundary condi- 
tions. ModelE2-R, like the multi-model PlioMIP ensemble, 
shows that the South Atlantic and southern Indian Ocean 
regions are somewhat warmer than the corresponding lati- 
tudes in the South Pacific. This feature tends to be robust, 
and is also borne out by the PRISM3D reconstructions. The 
results with GM.CORR appear to be an improvement com- 
pared to GM_UNCOR, which had broad warming throughout 
the southern high latitudes. However, ocean core data is very 
limited at many longitudes in the Southern Ocean and the 
meridional character of the Pliocene Southern Hemisphere 
sea surface temperatures is not well understood. 

Tropical temperature change, as mentioned previously, has 
proven to be a dilemma for many model-data comparisons 
for past warm periods in Earth history. One of the overar- 
ching themes from proxy studies of pre-Pleistocene warm 
climates is that tropical temperatures have been relatively 
stable compared to higher latitudes, resulting in decreased 


meridional temperature gradients (e.g., Zachos et al., 1994; 
Galfetti et al., 2007; Pearson et al., 2007; Dowsett et al., 
1996, 2010b). In contrast, climate model simulations of these 
past time periods show that when the forcing is either a well- 
mixed greenhouse gas or the result of changes to solar inso- 
lation, tropical temperatures respond measurably (e.g., Rind 
and Chandler, 1991; Chandler et al., 1994; Sloan et al., 1996; 
Haywood et al., 2000; Jiang et al., 2005). 

The mid-Pliocene, with its relatively lower levels of at- 
mospheric CO 2 (compared to earlier warm periods in geo- 
logic history), has minimal additional forcing in the tropics, 
thus, tropical SSTs and surface air temperatures warm only 
about 1 °C relative to the preindustrial. Tropical temperatures 
warmed another a degree in the Atlantic Ocean when the 
ocean mixing was corrected, but that was related to the fact 
that the uncorrected model received substantial cool surface 
flow from the unrealistic cold anomaly in the North Atlantic. 
Otherwise there was little change in response in the tropics 
between GM_CORR and GMTJNCORR. 

Results from the PRISM analyses show that some tropi- 
cal regions may in fact have experienced warming on the or- 
der of 1-2 °C, but PRISM also includes a number of sites 
where no temperature increases are discernable. In each 
ocean basin, the sites that show little or no temperature 
change are in the central or western warm pools, whereas 
those sites with minor warming are generally in the eastern- 
most portions of the basins. All of the ModelE2-R simula- 
tion results show generally uniform east-west temperature in- 
creases across all ocean basins and, where minor variations 
exist, there is no consistent east-west trend. This permanent 
El Nino or “El Padre” character of the tropical SSTs is con- 
sistent with some Pliocene data and simulations (Shukla et 
al., 2009 and references therein), but the issue is far from 
resolved with both data and modelling studies suggesting 
ENSO variability was still present in the Pliocene (Scroxton 
et al., 2011; Watanabe et al., 2011). 

3.5 Arctic sea ice 

The impact of the Pliocene warming on Arctic sea ice is of 
special interest because of the continuing degradation of the 
present-day polar cryosphere related to greenhouse warm- 
ing. Coupled model experiments calculate rather than spec- 
ify the sea ice distribution, and although the initial conditions 
for these simulations include reduced ice, it is not clear that 
models can maintain such distributions with the 405 ppmv 
CO 2 levels used in the PlioMIP experiments. In the case 
of the GISS ModelE2-R experiments, however, the reduc- 
tion is at least partially sustained: sea ice is reduced in both 
the GM.CORR and GM_UNCOR versions (Fig. 4 and Ta- 
ble 3). The corrected ocean mixing scheme shows sea ice 
is reduced by 40 % in the Arctic, which is nearly twice as 
before the correction was applied. (Note that in the South- 
ern Hemisphere (maps not shown), the Pliocene ocean ice is 
also less compared to the preindustrial, but the ocean mixing 
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SEA ICE COVERAGE 




Fig. 4. Extent of sea ice cover from the two preindustrial simulations (top) and the mid-Pliocene simulations (bottom) for both the 
GM_UNCOR and GM_CORR versions of the model. Sea ice in the corrected Pliocene run is reduced by 40 % in the Arctic Ocean and 
retreats on all margins, not just in the North Atlantic where SSTs warm the most. The reduction of ice was less before the ocean mixing 
scheme correction was applied, but the geographic character of ice reduction is similar. The central Arctic ice cap in the Pliocene no longer 
maintains a large geographic area where ice cover is above 90 %. Ice thickness is greatly reduced as well, indicative of a loss of multi-year 
ice and suggesting that the region would be highly vulnerable to even moderate further increases in temperature. 


correction actually reduces the amount of ice lost). Regard- 
less, the geographic character of the reduction is very sim- 
ilar between both corrected and uncorrected versions of the 
model in the Arctic. Reductions of sea ice are greatest in the 
proximity of the GIN Sea warming, but decreased sea ice 
extent and thickness are pronounced everywhere around the 
Arctic ice cap, especially at the margins of the ice cap. 

A further examination of the vertical temperature profile 
in the Atlantic and Arctic sector (Fig. 5) also shows that 
the remaining Arctic sea ice may be extremely vulnerable to 
even small amounts of additional warming. The Arctic Ocean 
warms considerably at depth and the sea ice cap thins dramat- 
ically as its areal coverage is reduced. No doubt these results 
are highly sensitive to the ocean-ice parameterisation in the 


climate model. We explore that issue, along with the seasonal 
cycle of polar sea ice, in greater detail in a forthcoming pa- 
per where we will discuss further the impact of the corrected 
ocean mixing scheme on the deep ocean. 

3.6 Atlantic meridional overturning circulation 

One of the most significant differences of the Pliocene 
GM_CORR simulations, compared with those of the uncor- 
rected model, is the characteristic of the meridional over- 
turning in the Atlantic Ocean. In GM_UNCOR the Atlantic 
Meridional Overturning Circulation (AMOC) collapsed and 
did not recover, something that was expected to be related 
to problems with the ocean mixing scheme. Although we 
hesitate to state that this is a clear improvement (little direct 
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Atlantic Vertical Temperature Profile 


Preindustrial Control -1.95 28.85 Pliocene -1.24 29.87 
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Fig. 5. The vertical profile of zonal average ocean temperatures 
in the Atlantic Ocean basin for the corrected (GM_CORR) Prein- 
dustrial control run (left) and the Pliocene simulation (right). Deep 
and bottom water temperatures warm, associated with reductions in 
Antarctic Bottom Water production (not shown) and increased pro- 
duction of North Atlantic Deep Water (see Fig. 6 and text). Warming 
at depth in the Arctic Ocean is also evident, leaving sea ice at the 
surface even more vulnerable to melting than the substantial 40 % 
ice cover reduction might suggest. 

evidence from observations), it seems likely that the col- 
lapsed AMOC in the previous simulation was erroneous. It 
is also uncertain how the 25 % increase in the AMOC con- 
tributes to the warmer SSTs since the associated increase in 
ocean heat transports is only 4 % (Zhang et al., 2013). 

The Atlantic mass stream function for GM_CORR is 
shown in Fig. 6 for both the preindustrial control and the 
Pliocene simulation. The peak overturning increases from 
14.35 to 17.89 Sverdrups, a 25 % strengthening of the cir- 
culation, and the geographic location of the peak shifts pole- 
ward by 5-10° latitude. Also apparent is the formation of an 
overturning cell at about 70° N, which appears to derive deep 
water production from either the anomalously warm Iceland 
or Norwegian Seas (or both). Taken together with the ver- 
tical temperature profile shown in Fig. 5, this newly devel- 
oped overturning cell may be the origin of warming ocean 
waters at depth in the Arctic Ocean. Further sensitivity ex- 
periments will be required to examine to what extent these 
features are driven by (1) the increased CO 2 , (2) the reduc- 
tion of the Greenland Ice Sheet, (3) changes in freshwater 
flux associated with altered river outflow, and (4) what affect 
the 25-m sea level rise may have had on flow along coastlines 
or through sub grid-scale straits. 


Atlantic Mass Stream Function 



(Sv) 
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Fig. 6. Zonally averaged mass stream function in the Atlantic 
Ocean basin for the corrected (GM_CORR) Preindustrial (left) and 
Pliocene simulations (right). Peak overturning in the North Atlantic 
of the Pliocene strengthens by over 25 % compared to the preindus- 
trial run, including an intrusion of a plume of water, with origins in 
the North Atlantic, to depths of over 1000 m in the Arctic Ocean. 
This increased overturning in the Pliocene has been hypothesised 
for years based on surface and deep water temperature proxies and 
estimates of NADW penetration relative to water masses from other 
basins. However, coupled climate models have rarely found this fea- 
ture in Pliocene experiments. 


4 Hydrology 

In most warm climate scenarios, the hydrological cycle 
strengthens as more water is evaporated from the oceans and 
the water holding capacity of the atmosphere increases. Sur- 
face heating destabilises air masses, increasing convective 
events leading to increases in precipitation rate, punctuated 
by longer periods of drying over land. Table 3 shows some of 
these effects on global hydrological variables as atmospheric 
water vapour, precipitation, evaporation and moist convec- 
tion all increase in the Pliocene. Total cloud cover actually 
decreases slightly, a characteristic of some warm climate 
simulations (Wang and Lau, 2006; Haywood et al., 2009; 
Bender, 2011) and supported by recent observations (Liu et 
al., 2009; Tang et al., 2012). The ocean mixing scheme cor- 
rection has, somewhat surprisingly, little effect on the hydro- 
logical cycle globally averaged. As mentioned, the Pliocene 
precipitation and evaporation rates increase globally in both 
GMJJNCOR and GM_CORR relative to the preindustrial. 
However, precipitation over continents in GM.CORR actu- 
ally increases more than evaporation (see Fig. 7), leading to 
a net moisture gain of approximately 1 5 % over land. At the 
same time, precipitation and evaporation anomalies over land 
are highly variable and local effects dominate the moisture 
balance. The vast majority of continental regions show that 
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increased rainfall is offset by increased evaporation (com- 
pare Fig. 7b to c); similarly, decreasing precipitation over 
land is commonly accompanied by lower evaporation rates 
in the same region. Despite the distinct differences when 
the Pliocene is compared with the preindustrial, the similar- 
ity between both the uncorrected and corrected versions of 
ModelE2-R are striking. The only major differences that are 
apparent in the precipitation anomalies occur over the North 
Atlantic, where the uncorrected ocean mixing error leads to 
cooling temperatures (even compared to the preindustrial) 
and a dampened hydrological cycle. The precipitation and 
evaporation fields for GM.CORR are fully consistent with 
the SST anomalies shown in Fig. 3. 

Over ocean regions, however, this is not necessarily true. 
In particular, the subtropical oceans experience intensified 
drying, as regions of increased evaporation are commonly 
affected by declining precipitation (shown for GM_CORR 
only). This effect is seen in the subtropics of the Pacific, At- 
lantic and Indian Oceans, as well as in both hemispheres. 
The implications for salinity distributions and the meridional 
overturning circulation may be significant. We are in the pro- 
cess of examining the impacts, and will discuss the findings 
in a future paper as further sensitivity experiments are com- 
pleted. 


5 Discussion of global feedbacks 

The primary drivers behind the climate changes simulated 
in the ModelE2-R Pliocene simulations are threefold: (1) al- 
tered boundary conditions, including minor modifications 
to continental geography, orography and coastlines, (2) re- 
duced Greenland and West Antarctic ice sheets, and the pole- 
ward extension of the Southern Ocean over West Antarc- 
tica; and (3) an increased atmospheric carbon dioxide level 
(405 ppmv). However, as with all warm climate experiments, 
the amplifying effects of feedbacks are fundamental to the 
evolution of the simulation and to the ultimate equilibrium 
state of the climate. 

With global average temperatures increasing in the mid- 
Pliocene simulations more than 2°C compared to the 
pre industrial control runs and with the Northern Hemi- 
sphere warming more than the Southern Hemisphere (2.4 °C 
vs. 2.1 °C in GM.CORR) we look to changes in albedo and 
water vapour to explain much of the difference. All pro- 
vide important positive feedbacks in the warmer Pliocene cli- 
mates. The reduced planetary albedo in the Pliocene simula- 
tions is consistent across both Northern and Southern Hemi- 
spheres and is partly related to cloud cover reductions of a 
couple percent. Of course the ground albedo change is even 
more substantial (—20% NH, —13 % SH), as snow and ice 
cover decrease dramatically in the warmer climate. The spec- 
ified reduction of ice sheets imposes a decrease in ground 
albedo, although it is not as large an effect as the albedo 
change caused by the loss of the snow and ice cover, which 


PRECIPITATION 

Pliocene - Preindustrial (GM_UNCOR) 



Pliocene - Preindustrial (GM_CORR) 



EVAPORATION 

Pliocene - Preindustrial (GM_CORR) 



Fig. 7. Precipitation anomalies for both mid-Pliocene simula- 
tions minus their Preindustrial control runs (7a-GM_UNCOR; 7b- 
GM_CORR). The patterns are largely similar except in the North 
and tropical Atlantic where the ocean mixing correction substan- 
tially altered SSTs. An intensified hydrological cycle is evident 
from global increases in both rainfall and evaporation (7c-shown for 
GM.CORR only). Note that the colourbar is flipped for the evapo- 
ration map so that green shades on all maps indicate a tendency to- 
ward wetter conditions while brown indicates drier conditions. Most 
regional changes in precipitation over land are balanced by evapo- 
ration, at least in sign if not in magnitude. In contrast, there are 
significant regions of the ocean, especially in the subtropics, where 
the hydrological anomaly in both precipitation and evaporation am- 
plifies the other. The impact on ocean surface salinity needs to be 
explored with further experiments and analyses. 
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is calculated by the GCM. Snow depth, while decreasing in 
the Northern Hemisphere, actually increases on average in 
the Southern Hemisphere in the Pliocene simulations. This 
less obvious result is due to the great reduction in circum- 
Antarctic sea ice, which leads to increased evaporation across 
warm waters and large increases in snow accumulation over 
the margins of East Antarctica. Although the West Antarc- 
tica Ice Sheet has been removed, and West Antarctica it- 
self is mostly submerged, sea ice formation across that high 
latitude ocean still forms a platform that allows snow accu- 
mulation for much of the year. This begs the question as to 
whether or not the GISS ModelE2-R would actually support 
the re-initiation or maintenance of continental glaciation on 
any portion of West Antarctica. 

As mentioned above, the other factor in planetary albedo 
change is the role of cloud cover. Total cloud cover decreases 
in the Pliocene, but this is not necessarily an indication of the 
sign of the feedback, since the vertical distribution of clouds 
can alter their impact on albedo. However, the vertical pro- 
file of cloud changes shown in Fig. 8 strongly suggests that 
clouds enhance the Pliocene warming, since they decrease on 
average at every level in the lower and middle troposphere, 
while increasing slightly near the tropopause. Moist convec- 
tive cloud cover is the only cloud type to increase in the 
Pliocene simulations, consistent with the strengthened hy- 
drologic cycle. 

Beyond the albedo changes, the atmospheric water vapour 
increases by about 15 %, providing the typical amplification 
of the greenhouse effect and another positive feedback that 
maintains a warmer Pliocene climate. The increased surface 
area covered by water in the Pliocene plays a role in this in- 
creased atmospheric water vapour, but the negative moisture 
balance across the subtropical ocean expanses of the South- 
ern Hemisphere is a key factor as well (Fig. 7). Intensification 
of the hydrological cycle results from increases in both rain- 
fall and evaporation rates (+6 %), but annual average precip- 
itation changes over most continental regions are balanced 
by similarly altered regional patterns of evaporation (at least 
in sign if not entirely in magnitude). In contrast, there are 
significant geographic expanses of the subtropics where both 
precipitation and evaporation change amplifies the effect of 
the other. The impact on ocean surface salinity may be im- 
portant to the long-term equilibrium state of Pliocene ocean 
circulation; we are in the process of conducting longer in- 
tegrations (beyond these 1000 yr simulations) to explore the 
issue. 


6 Conclusions 

Given its intriguing potential as a near-analogue to future 
climate, the mid-Pliocene Warm Period has often been ex- 
plored by paleoclimatologists and modellers. The interval 
between 3. 3-3.0 million yr ago may truly be the most re- 
cent global warming that is anywhere close in magnitude to 


Vertical Cloud Cover Difference 



Fig. 8. The zonally averaged vertical profile of total cloud cover 
from the corrected (GM_CORR) Pliocene and Preindustrial simu- 
lations (left plot) and the difference between the two (right plot). 
Note the separate scales on the x-axes. Changes in cloud cover at 
all heights are small, but they are consistent. Total cloud cover de- 
creases throughout most of the troposphere and acts as a positive 
feedback to warming in the Pliocene simulation. 


what the Earth faces in its near future. Unfortunately, exper- 
iments with many of the coupled ocean-atmosphere models 
that IPCC relies on for future climate projections have not 
been fully successful at reproducing the middle Pliocene dis- 
tribution of warm sea surface temperatures, particularly the 
strong warming in the far North Atlantic Ocean (Haywood 
et al., 2013; Dowsett et al., 2012). Numerous proxy stud- 
ies imply that higher CO 2 levels are not the sole answer to 
Pliocene warming, as prior assertions of stable tropical SSTs 
have been corroborated with each new low-latitude ocean 
core and each new proxy method. 

In this paper, we show preliminary results of the most 
recent simulations of the mid-Pliocene Warm Period using 
the latest version (AR5/CMIP5) of the GISS Earth System 
Model, called ModelE2-R. We discuss two versions of the 
Pliocene and Preindustrial simulations because one set of 
simulations includes a post-CMIP5 correction to the model’s 
Gent-McWilliams ocean mixing scheme that has a substan- 
tial impact on the results - and offers a substantial im- 
provement, correcting some serious problems with the GISS 
ModelE2-R ocean. Both Pliocene simulations represent the 
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NASA/GISS contribution to the Pliocene Model Intercom- 
parison Project (PlioMIP, Experiment 2) and are used in 
Phase 1 of the PlioMIP multi-model ensemble studies. Go- 
ing forward, we will adopt only the model version that in- 
corporates the corrected mixing scheme (GM_CORR in this 
paper). Many of our results from the Pliocene GM_CORR 
simulation fall squarely within the range defined by other 
coupled models in the PlioMIP ensemble, but we empha- 
sise here some features that would be outliers in the en- 
semble ranges. The most prominent is the simulation of a 
large region of warming in the North Atlantic and Greenland- 
Iceland-Norwegian Sea (large areas covered by SST anoma- 
lies greater than +4°C, with distinct peaks of more than 
9°C). This is by far the most accurate portrayal of this key 
geographic region by any version of the NASA/GISS fam- 
ily of models to date and comparisons with the Pliocene 
GMJJNCOR simulations shows that the ocean mixing 
scheme correction is the key factor in the improved results 
- at least for ModelE2-R. There are still model-data differ- 
ences to be addressed, and these preliminary results require 
further analysis in order to fully understand the physical pro- 
cesses involved in the warming of the North Atlantic and in 
the changes to ocean circulation. However, we believe that 
continued development of key physical routines in the GISS 
atmospheric model, along with higher resolution and recent 
corrections made to the Gent-McWilliams mixing parame- 
terisation in the Russell ocean model, have led to an Earth 
System Model that will produce more accurate projections 
of both past and future climate. 
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